3.8.84 \(\int \frac {1}{x^2 (a+b x^4) (c+d x^4)} \, dx\) [784]

3.8.84.1 Optimal result
3.8.84.2 Mathematica [A] (verified)
3.8.84.3 Rubi [A] (verified)
3.8.84.4 Maple [A] (verified)
3.8.84.5 Fricas [C] (verification not implemented)
3.8.84.6 Sympy [F(-1)]
3.8.84.7 Maxima [A] (verification not implemented)
3.8.84.8 Giac [A] (verification not implemented)
3.8.84.9 Mupad [B] (verification not implemented)

3.8.84.1 Optimal result

Integrand size = 22, antiderivative size = 460 \[ \int \frac {1}{x^2 \left (a+b x^4\right ) \left (c+d x^4\right )} \, dx=-\frac {1}{a c x}+\frac {b^{5/4} \arctan \left (1-\frac {\sqrt {2} \sqrt [4]{b} x}{\sqrt [4]{a}}\right )}{2 \sqrt {2} a^{5/4} (b c-a d)}-\frac {b^{5/4} \arctan \left (1+\frac {\sqrt {2} \sqrt [4]{b} x}{\sqrt [4]{a}}\right )}{2 \sqrt {2} a^{5/4} (b c-a d)}-\frac {d^{5/4} \arctan \left (1-\frac {\sqrt {2} \sqrt [4]{d} x}{\sqrt [4]{c}}\right )}{2 \sqrt {2} c^{5/4} (b c-a d)}+\frac {d^{5/4} \arctan \left (1+\frac {\sqrt {2} \sqrt [4]{d} x}{\sqrt [4]{c}}\right )}{2 \sqrt {2} c^{5/4} (b c-a d)}-\frac {b^{5/4} \log \left (\sqrt {a}-\sqrt {2} \sqrt [4]{a} \sqrt [4]{b} x+\sqrt {b} x^2\right )}{4 \sqrt {2} a^{5/4} (b c-a d)}+\frac {b^{5/4} \log \left (\sqrt {a}+\sqrt {2} \sqrt [4]{a} \sqrt [4]{b} x+\sqrt {b} x^2\right )}{4 \sqrt {2} a^{5/4} (b c-a d)}+\frac {d^{5/4} \log \left (\sqrt {c}-\sqrt {2} \sqrt [4]{c} \sqrt [4]{d} x+\sqrt {d} x^2\right )}{4 \sqrt {2} c^{5/4} (b c-a d)}-\frac {d^{5/4} \log \left (\sqrt {c}+\sqrt {2} \sqrt [4]{c} \sqrt [4]{d} x+\sqrt {d} x^2\right )}{4 \sqrt {2} c^{5/4} (b c-a d)} \]

output
-1/a/c/x-1/4*b^(5/4)*arctan(-1+b^(1/4)*x*2^(1/2)/a^(1/4))/a^(5/4)/(-a*d+b* 
c)*2^(1/2)-1/4*b^(5/4)*arctan(1+b^(1/4)*x*2^(1/2)/a^(1/4))/a^(5/4)/(-a*d+b 
*c)*2^(1/2)+1/4*d^(5/4)*arctan(-1+d^(1/4)*x*2^(1/2)/c^(1/4))/c^(5/4)/(-a*d 
+b*c)*2^(1/2)+1/4*d^(5/4)*arctan(1+d^(1/4)*x*2^(1/2)/c^(1/4))/c^(5/4)/(-a* 
d+b*c)*2^(1/2)-1/8*b^(5/4)*ln(-a^(1/4)*b^(1/4)*x*2^(1/2)+a^(1/2)+x^2*b^(1/ 
2))/a^(5/4)/(-a*d+b*c)*2^(1/2)+1/8*b^(5/4)*ln(a^(1/4)*b^(1/4)*x*2^(1/2)+a^ 
(1/2)+x^2*b^(1/2))/a^(5/4)/(-a*d+b*c)*2^(1/2)+1/8*d^(5/4)*ln(-c^(1/4)*d^(1 
/4)*x*2^(1/2)+c^(1/2)+x^2*d^(1/2))/c^(5/4)/(-a*d+b*c)*2^(1/2)-1/8*d^(5/4)* 
ln(c^(1/4)*d^(1/4)*x*2^(1/2)+c^(1/2)+x^2*d^(1/2))/c^(5/4)/(-a*d+b*c)*2^(1/ 
2)
 
3.8.84.2 Mathematica [A] (verified)

Time = 0.21 (sec) , antiderivative size = 385, normalized size of antiderivative = 0.84 \[ \int \frac {1}{x^2 \left (a+b x^4\right ) \left (c+d x^4\right )} \, dx=\frac {\frac {8 b}{a}-\frac {8 d}{c}-\frac {2 \sqrt {2} b^{5/4} x \arctan \left (1-\frac {\sqrt {2} \sqrt [4]{b} x}{\sqrt [4]{a}}\right )}{a^{5/4}}+\frac {2 \sqrt {2} b^{5/4} x \arctan \left (1+\frac {\sqrt {2} \sqrt [4]{b} x}{\sqrt [4]{a}}\right )}{a^{5/4}}+\frac {2 \sqrt {2} d^{5/4} x \arctan \left (1-\frac {\sqrt {2} \sqrt [4]{d} x}{\sqrt [4]{c}}\right )}{c^{5/4}}-\frac {2 \sqrt {2} d^{5/4} x \arctan \left (1+\frac {\sqrt {2} \sqrt [4]{d} x}{\sqrt [4]{c}}\right )}{c^{5/4}}+\frac {\sqrt {2} b^{5/4} x \log \left (\sqrt {a}-\sqrt {2} \sqrt [4]{a} \sqrt [4]{b} x+\sqrt {b} x^2\right )}{a^{5/4}}-\frac {\sqrt {2} b^{5/4} x \log \left (\sqrt {a}+\sqrt {2} \sqrt [4]{a} \sqrt [4]{b} x+\sqrt {b} x^2\right )}{a^{5/4}}-\frac {\sqrt {2} d^{5/4} x \log \left (\sqrt {c}-\sqrt {2} \sqrt [4]{c} \sqrt [4]{d} x+\sqrt {d} x^2\right )}{c^{5/4}}+\frac {\sqrt {2} d^{5/4} x \log \left (\sqrt {c}+\sqrt {2} \sqrt [4]{c} \sqrt [4]{d} x+\sqrt {d} x^2\right )}{c^{5/4}}}{-8 b c x+8 a d x} \]

input
Integrate[1/(x^2*(a + b*x^4)*(c + d*x^4)),x]
 
output
((8*b)/a - (8*d)/c - (2*Sqrt[2]*b^(5/4)*x*ArcTan[1 - (Sqrt[2]*b^(1/4)*x)/a 
^(1/4)])/a^(5/4) + (2*Sqrt[2]*b^(5/4)*x*ArcTan[1 + (Sqrt[2]*b^(1/4)*x)/a^( 
1/4)])/a^(5/4) + (2*Sqrt[2]*d^(5/4)*x*ArcTan[1 - (Sqrt[2]*d^(1/4)*x)/c^(1/ 
4)])/c^(5/4) - (2*Sqrt[2]*d^(5/4)*x*ArcTan[1 + (Sqrt[2]*d^(1/4)*x)/c^(1/4) 
])/c^(5/4) + (Sqrt[2]*b^(5/4)*x*Log[Sqrt[a] - Sqrt[2]*a^(1/4)*b^(1/4)*x + 
Sqrt[b]*x^2])/a^(5/4) - (Sqrt[2]*b^(5/4)*x*Log[Sqrt[a] + Sqrt[2]*a^(1/4)*b 
^(1/4)*x + Sqrt[b]*x^2])/a^(5/4) - (Sqrt[2]*d^(5/4)*x*Log[Sqrt[c] - Sqrt[2 
]*c^(1/4)*d^(1/4)*x + Sqrt[d]*x^2])/c^(5/4) + (Sqrt[2]*d^(5/4)*x*Log[Sqrt[ 
c] + Sqrt[2]*c^(1/4)*d^(1/4)*x + Sqrt[d]*x^2])/c^(5/4))/(-8*b*c*x + 8*a*d* 
x)
 
3.8.84.3 Rubi [A] (verified)

Time = 0.64 (sec) , antiderivative size = 477, normalized size of antiderivative = 1.04, number of steps used = 4, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.182, Rules used = {980, 25, 1054, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {1}{x^2 \left (a+b x^4\right ) \left (c+d x^4\right )} \, dx\)

\(\Big \downarrow \) 980

\(\displaystyle \frac {\int -\frac {x^2 \left (b d x^4+b c+a d\right )}{\left (b x^4+a\right ) \left (d x^4+c\right )}dx}{a c}-\frac {1}{a c x}\)

\(\Big \downarrow \) 25

\(\displaystyle -\frac {\int \frac {x^2 \left (b d x^4+b c+a d\right )}{\left (b x^4+a\right ) \left (d x^4+c\right )}dx}{a c}-\frac {1}{a c x}\)

\(\Big \downarrow \) 1054

\(\displaystyle -\frac {\int \left (\frac {b^2 c x^2}{(b c-a d) \left (b x^4+a\right )}+\frac {a d^2 x^2}{(a d-b c) \left (d x^4+c\right )}\right )dx}{a c}-\frac {1}{a c x}\)

\(\Big \downarrow \) 2009

\(\displaystyle -\frac {-\frac {b^{5/4} c \arctan \left (1-\frac {\sqrt {2} \sqrt [4]{b} x}{\sqrt [4]{a}}\right )}{2 \sqrt {2} \sqrt [4]{a} (b c-a d)}+\frac {b^{5/4} c \arctan \left (\frac {\sqrt {2} \sqrt [4]{b} x}{\sqrt [4]{a}}+1\right )}{2 \sqrt {2} \sqrt [4]{a} (b c-a d)}+\frac {a d^{5/4} \arctan \left (1-\frac {\sqrt {2} \sqrt [4]{d} x}{\sqrt [4]{c}}\right )}{2 \sqrt {2} \sqrt [4]{c} (b c-a d)}-\frac {a d^{5/4} \arctan \left (\frac {\sqrt {2} \sqrt [4]{d} x}{\sqrt [4]{c}}+1\right )}{2 \sqrt {2} \sqrt [4]{c} (b c-a d)}+\frac {b^{5/4} c \log \left (-\sqrt {2} \sqrt [4]{a} \sqrt [4]{b} x+\sqrt {a}+\sqrt {b} x^2\right )}{4 \sqrt {2} \sqrt [4]{a} (b c-a d)}-\frac {b^{5/4} c \log \left (\sqrt {2} \sqrt [4]{a} \sqrt [4]{b} x+\sqrt {a}+\sqrt {b} x^2\right )}{4 \sqrt {2} \sqrt [4]{a} (b c-a d)}-\frac {a d^{5/4} \log \left (-\sqrt {2} \sqrt [4]{c} \sqrt [4]{d} x+\sqrt {c}+\sqrt {d} x^2\right )}{4 \sqrt {2} \sqrt [4]{c} (b c-a d)}+\frac {a d^{5/4} \log \left (\sqrt {2} \sqrt [4]{c} \sqrt [4]{d} x+\sqrt {c}+\sqrt {d} x^2\right )}{4 \sqrt {2} \sqrt [4]{c} (b c-a d)}}{a c}-\frac {1}{a c x}\)

input
Int[1/(x^2*(a + b*x^4)*(c + d*x^4)),x]
 
output
-(1/(a*c*x)) - (-1/2*(b^(5/4)*c*ArcTan[1 - (Sqrt[2]*b^(1/4)*x)/a^(1/4)])/( 
Sqrt[2]*a^(1/4)*(b*c - a*d)) + (b^(5/4)*c*ArcTan[1 + (Sqrt[2]*b^(1/4)*x)/a 
^(1/4)])/(2*Sqrt[2]*a^(1/4)*(b*c - a*d)) + (a*d^(5/4)*ArcTan[1 - (Sqrt[2]* 
d^(1/4)*x)/c^(1/4)])/(2*Sqrt[2]*c^(1/4)*(b*c - a*d)) - (a*d^(5/4)*ArcTan[1 
 + (Sqrt[2]*d^(1/4)*x)/c^(1/4)])/(2*Sqrt[2]*c^(1/4)*(b*c - a*d)) + (b^(5/4 
)*c*Log[Sqrt[a] - Sqrt[2]*a^(1/4)*b^(1/4)*x + Sqrt[b]*x^2])/(4*Sqrt[2]*a^( 
1/4)*(b*c - a*d)) - (b^(5/4)*c*Log[Sqrt[a] + Sqrt[2]*a^(1/4)*b^(1/4)*x + S 
qrt[b]*x^2])/(4*Sqrt[2]*a^(1/4)*(b*c - a*d)) - (a*d^(5/4)*Log[Sqrt[c] - Sq 
rt[2]*c^(1/4)*d^(1/4)*x + Sqrt[d]*x^2])/(4*Sqrt[2]*c^(1/4)*(b*c - a*d)) + 
(a*d^(5/4)*Log[Sqrt[c] + Sqrt[2]*c^(1/4)*d^(1/4)*x + Sqrt[d]*x^2])/(4*Sqrt 
[2]*c^(1/4)*(b*c - a*d)))/(a*c)
 

3.8.84.3.1 Defintions of rubi rules used

rule 25
Int[-(Fx_), x_Symbol] :> Simp[Identity[-1]   Int[Fx, x], x]
 

rule 980
Int[((e_.)*(x_))^(m_)*((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_) 
)^(q_), x_Symbol] :> Simp[(e*x)^(m + 1)*(a + b*x^n)^(p + 1)*((c + d*x^n)^(q 
 + 1)/(a*c*e*(m + 1))), x] - Simp[1/(a*c*e^n*(m + 1))   Int[(e*x)^(m + n)*( 
a + b*x^n)^p*(c + d*x^n)^q*Simp[(b*c + a*d)*(m + n + 1) + n*(b*c*p + a*d*q) 
 + b*d*(m + n*(p + q + 2) + 1)*x^n, x], x], x] /; FreeQ[{a, b, c, d, e, p, 
q}, x] && NeQ[b*c - a*d, 0] && IGtQ[n, 0] && LtQ[m, -1] && IntBinomialQ[a, 
b, c, d, e, m, n, p, q, x]
 

rule 1054
Int[(((g_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_)*((e_) + (f_.)*(x_)^(n 
_)))/((c_) + (d_.)*(x_)^(n_)), x_Symbol] :> Int[ExpandIntegrand[(g*x)^m*(a 
+ b*x^n)^p*((e + f*x^n)/(c + d*x^n)), x], x] /; FreeQ[{a, b, c, d, e, f, g, 
 m, p}, x] && IGtQ[n, 0]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 
3.8.84.4 Maple [A] (verified)

Time = 4.68 (sec) , antiderivative size = 237, normalized size of antiderivative = 0.52

method result size
default \(\frac {b \sqrt {2}\, \left (\ln \left (\frac {x^{2}-\left (\frac {a}{b}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {a}{b}}}{x^{2}+\left (\frac {a}{b}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {a}{b}}}\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {a}{b}\right )^{\frac {1}{4}}}+1\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {a}{b}\right )^{\frac {1}{4}}}-1\right )\right )}{8 \left (a d -b c \right ) a \left (\frac {a}{b}\right )^{\frac {1}{4}}}-\frac {1}{a c x}-\frac {d \sqrt {2}\, \left (\ln \left (\frac {x^{2}-\left (\frac {c}{d}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {c}{d}}}{x^{2}+\left (\frac {c}{d}\right )^{\frac {1}{4}} x \sqrt {2}+\sqrt {\frac {c}{d}}}\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {c}{d}\right )^{\frac {1}{4}}}+1\right )+2 \arctan \left (\frac {\sqrt {2}\, x}{\left (\frac {c}{d}\right )^{\frac {1}{4}}}-1\right )\right )}{8 \left (a d -b c \right ) c \left (\frac {c}{d}\right )^{\frac {1}{4}}}\) \(237\)
risch \(\text {Expression too large to display}\) \(1027\)

input
int(1/x^2/(b*x^4+a)/(d*x^4+c),x,method=_RETURNVERBOSE)
 
output
1/8*b/(a*d-b*c)/a/(a/b)^(1/4)*2^(1/2)*(ln((x^2-(a/b)^(1/4)*x*2^(1/2)+(a/b) 
^(1/2))/(x^2+(a/b)^(1/4)*x*2^(1/2)+(a/b)^(1/2)))+2*arctan(2^(1/2)/(a/b)^(1 
/4)*x+1)+2*arctan(2^(1/2)/(a/b)^(1/4)*x-1))-1/a/c/x-1/8*d/(a*d-b*c)/c/(c/d 
)^(1/4)*2^(1/2)*(ln((x^2-(c/d)^(1/4)*x*2^(1/2)+(c/d)^(1/2))/(x^2+(c/d)^(1/ 
4)*x*2^(1/2)+(c/d)^(1/2)))+2*arctan(2^(1/2)/(c/d)^(1/4)*x+1)+2*arctan(2^(1 
/2)/(c/d)^(1/4)*x-1))
 
3.8.84.5 Fricas [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.50 (sec) , antiderivative size = 1461, normalized size of antiderivative = 3.18 \[ \int \frac {1}{x^2 \left (a+b x^4\right ) \left (c+d x^4\right )} \, dx=\text {Too large to display} \]

input
integrate(1/x^2/(b*x^4+a)/(d*x^4+c),x, algorithm="fricas")
 
output
-1/4*((-b^5/(a^5*b^4*c^4 - 4*a^6*b^3*c^3*d + 6*a^7*b^2*c^2*d^2 - 4*a^8*b*c 
*d^3 + a^9*d^4))^(1/4)*a*c*x*log(b^4*x + (a^4*b^3*c^3 - 3*a^5*b^2*c^2*d + 
3*a^6*b*c*d^2 - a^7*d^3)*(-b^5/(a^5*b^4*c^4 - 4*a^6*b^3*c^3*d + 6*a^7*b^2* 
c^2*d^2 - 4*a^8*b*c*d^3 + a^9*d^4))^(3/4)) - (-b^5/(a^5*b^4*c^4 - 4*a^6*b^ 
3*c^3*d + 6*a^7*b^2*c^2*d^2 - 4*a^8*b*c*d^3 + a^9*d^4))^(1/4)*a*c*x*log(b^ 
4*x - (a^4*b^3*c^3 - 3*a^5*b^2*c^2*d + 3*a^6*b*c*d^2 - a^7*d^3)*(-b^5/(a^5 
*b^4*c^4 - 4*a^6*b^3*c^3*d + 6*a^7*b^2*c^2*d^2 - 4*a^8*b*c*d^3 + a^9*d^4)) 
^(3/4)) + I*(-b^5/(a^5*b^4*c^4 - 4*a^6*b^3*c^3*d + 6*a^7*b^2*c^2*d^2 - 4*a 
^8*b*c*d^3 + a^9*d^4))^(1/4)*a*c*x*log(b^4*x - (I*a^4*b^3*c^3 - 3*I*a^5*b^ 
2*c^2*d + 3*I*a^6*b*c*d^2 - I*a^7*d^3)*(-b^5/(a^5*b^4*c^4 - 4*a^6*b^3*c^3* 
d + 6*a^7*b^2*c^2*d^2 - 4*a^8*b*c*d^3 + a^9*d^4))^(3/4)) - I*(-b^5/(a^5*b^ 
4*c^4 - 4*a^6*b^3*c^3*d + 6*a^7*b^2*c^2*d^2 - 4*a^8*b*c*d^3 + a^9*d^4))^(1 
/4)*a*c*x*log(b^4*x - (-I*a^4*b^3*c^3 + 3*I*a^5*b^2*c^2*d - 3*I*a^6*b*c*d^ 
2 + I*a^7*d^3)*(-b^5/(a^5*b^4*c^4 - 4*a^6*b^3*c^3*d + 6*a^7*b^2*c^2*d^2 - 
4*a^8*b*c*d^3 + a^9*d^4))^(3/4)) - (-d^5/(b^4*c^9 - 4*a*b^3*c^8*d + 6*a^2* 
b^2*c^7*d^2 - 4*a^3*b*c^6*d^3 + a^4*c^5*d^4))^(1/4)*a*c*x*log(d^4*x + (b^3 
*c^7 - 3*a*b^2*c^6*d + 3*a^2*b*c^5*d^2 - a^3*c^4*d^3)*(-d^5/(b^4*c^9 - 4*a 
*b^3*c^8*d + 6*a^2*b^2*c^7*d^2 - 4*a^3*b*c^6*d^3 + a^4*c^5*d^4))^(3/4)) + 
(-d^5/(b^4*c^9 - 4*a*b^3*c^8*d + 6*a^2*b^2*c^7*d^2 - 4*a^3*b*c^6*d^3 + a^4 
*c^5*d^4))^(1/4)*a*c*x*log(d^4*x - (b^3*c^7 - 3*a*b^2*c^6*d + 3*a^2*b*c...
 
3.8.84.6 Sympy [F(-1)]

Timed out. \[ \int \frac {1}{x^2 \left (a+b x^4\right ) \left (c+d x^4\right )} \, dx=\text {Timed out} \]

input
integrate(1/x**2/(b*x**4+a)/(d*x**4+c),x)
 
output
Timed out
 
3.8.84.7 Maxima [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 384, normalized size of antiderivative = 0.83 \[ \int \frac {1}{x^2 \left (a+b x^4\right ) \left (c+d x^4\right )} \, dx=-\frac {b^{2} {\left (\frac {2 \, \sqrt {2} \arctan \left (\frac {\sqrt {2} {\left (2 \, \sqrt {b} x + \sqrt {2} a^{\frac {1}{4}} b^{\frac {1}{4}}\right )}}{2 \, \sqrt {\sqrt {a} \sqrt {b}}}\right )}{\sqrt {\sqrt {a} \sqrt {b}} \sqrt {b}} + \frac {2 \, \sqrt {2} \arctan \left (\frac {\sqrt {2} {\left (2 \, \sqrt {b} x - \sqrt {2} a^{\frac {1}{4}} b^{\frac {1}{4}}\right )}}{2 \, \sqrt {\sqrt {a} \sqrt {b}}}\right )}{\sqrt {\sqrt {a} \sqrt {b}} \sqrt {b}} - \frac {\sqrt {2} \log \left (\sqrt {b} x^{2} + \sqrt {2} a^{\frac {1}{4}} b^{\frac {1}{4}} x + \sqrt {a}\right )}{a^{\frac {1}{4}} b^{\frac {3}{4}}} + \frac {\sqrt {2} \log \left (\sqrt {b} x^{2} - \sqrt {2} a^{\frac {1}{4}} b^{\frac {1}{4}} x + \sqrt {a}\right )}{a^{\frac {1}{4}} b^{\frac {3}{4}}}\right )}}{8 \, {\left (a b c - a^{2} d\right )}} + \frac {d^{2} {\left (\frac {2 \, \sqrt {2} \arctan \left (\frac {\sqrt {2} {\left (2 \, \sqrt {d} x + \sqrt {2} c^{\frac {1}{4}} d^{\frac {1}{4}}\right )}}{2 \, \sqrt {\sqrt {c} \sqrt {d}}}\right )}{\sqrt {\sqrt {c} \sqrt {d}} \sqrt {d}} + \frac {2 \, \sqrt {2} \arctan \left (\frac {\sqrt {2} {\left (2 \, \sqrt {d} x - \sqrt {2} c^{\frac {1}{4}} d^{\frac {1}{4}}\right )}}{2 \, \sqrt {\sqrt {c} \sqrt {d}}}\right )}{\sqrt {\sqrt {c} \sqrt {d}} \sqrt {d}} - \frac {\sqrt {2} \log \left (\sqrt {d} x^{2} + \sqrt {2} c^{\frac {1}{4}} d^{\frac {1}{4}} x + \sqrt {c}\right )}{c^{\frac {1}{4}} d^{\frac {3}{4}}} + \frac {\sqrt {2} \log \left (\sqrt {d} x^{2} - \sqrt {2} c^{\frac {1}{4}} d^{\frac {1}{4}} x + \sqrt {c}\right )}{c^{\frac {1}{4}} d^{\frac {3}{4}}}\right )}}{8 \, {\left (b c^{2} - a c d\right )}} - \frac {1}{a c x} \]

input
integrate(1/x^2/(b*x^4+a)/(d*x^4+c),x, algorithm="maxima")
 
output
-1/8*b^2*(2*sqrt(2)*arctan(1/2*sqrt(2)*(2*sqrt(b)*x + sqrt(2)*a^(1/4)*b^(1 
/4))/sqrt(sqrt(a)*sqrt(b)))/(sqrt(sqrt(a)*sqrt(b))*sqrt(b)) + 2*sqrt(2)*ar 
ctan(1/2*sqrt(2)*(2*sqrt(b)*x - sqrt(2)*a^(1/4)*b^(1/4))/sqrt(sqrt(a)*sqrt 
(b)))/(sqrt(sqrt(a)*sqrt(b))*sqrt(b)) - sqrt(2)*log(sqrt(b)*x^2 + sqrt(2)* 
a^(1/4)*b^(1/4)*x + sqrt(a))/(a^(1/4)*b^(3/4)) + sqrt(2)*log(sqrt(b)*x^2 - 
 sqrt(2)*a^(1/4)*b^(1/4)*x + sqrt(a))/(a^(1/4)*b^(3/4)))/(a*b*c - a^2*d) + 
 1/8*d^2*(2*sqrt(2)*arctan(1/2*sqrt(2)*(2*sqrt(d)*x + sqrt(2)*c^(1/4)*d^(1 
/4))/sqrt(sqrt(c)*sqrt(d)))/(sqrt(sqrt(c)*sqrt(d))*sqrt(d)) + 2*sqrt(2)*ar 
ctan(1/2*sqrt(2)*(2*sqrt(d)*x - sqrt(2)*c^(1/4)*d^(1/4))/sqrt(sqrt(c)*sqrt 
(d)))/(sqrt(sqrt(c)*sqrt(d))*sqrt(d)) - sqrt(2)*log(sqrt(d)*x^2 + sqrt(2)* 
c^(1/4)*d^(1/4)*x + sqrt(c))/(c^(1/4)*d^(3/4)) + sqrt(2)*log(sqrt(d)*x^2 - 
 sqrt(2)*c^(1/4)*d^(1/4)*x + sqrt(c))/(c^(1/4)*d^(3/4)))/(b*c^2 - a*c*d) - 
 1/(a*c*x)
 
3.8.84.8 Giac [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 488, normalized size of antiderivative = 1.06 \[ \int \frac {1}{x^2 \left (a+b x^4\right ) \left (c+d x^4\right )} \, dx=-\frac {\left (a b^{3}\right )^{\frac {3}{4}} \arctan \left (\frac {\sqrt {2} {\left (2 \, x + \sqrt {2} \left (\frac {a}{b}\right )^{\frac {1}{4}}\right )}}{2 \, \left (\frac {a}{b}\right )^{\frac {1}{4}}}\right )}{2 \, {\left (\sqrt {2} a^{2} b^{2} c - \sqrt {2} a^{3} b d\right )}} - \frac {\left (a b^{3}\right )^{\frac {3}{4}} \arctan \left (\frac {\sqrt {2} {\left (2 \, x - \sqrt {2} \left (\frac {a}{b}\right )^{\frac {1}{4}}\right )}}{2 \, \left (\frac {a}{b}\right )^{\frac {1}{4}}}\right )}{2 \, {\left (\sqrt {2} a^{2} b^{2} c - \sqrt {2} a^{3} b d\right )}} + \frac {\left (c d^{3}\right )^{\frac {3}{4}} \arctan \left (\frac {\sqrt {2} {\left (2 \, x + \sqrt {2} \left (\frac {c}{d}\right )^{\frac {1}{4}}\right )}}{2 \, \left (\frac {c}{d}\right )^{\frac {1}{4}}}\right )}{2 \, {\left (\sqrt {2} b c^{3} d - \sqrt {2} a c^{2} d^{2}\right )}} + \frac {\left (c d^{3}\right )^{\frac {3}{4}} \arctan \left (\frac {\sqrt {2} {\left (2 \, x - \sqrt {2} \left (\frac {c}{d}\right )^{\frac {1}{4}}\right )}}{2 \, \left (\frac {c}{d}\right )^{\frac {1}{4}}}\right )}{2 \, {\left (\sqrt {2} b c^{3} d - \sqrt {2} a c^{2} d^{2}\right )}} + \frac {\left (a b^{3}\right )^{\frac {3}{4}} \log \left (x^{2} + \sqrt {2} x \left (\frac {a}{b}\right )^{\frac {1}{4}} + \sqrt {\frac {a}{b}}\right )}{4 \, {\left (\sqrt {2} a^{2} b^{2} c - \sqrt {2} a^{3} b d\right )}} - \frac {\left (a b^{3}\right )^{\frac {3}{4}} \log \left (x^{2} - \sqrt {2} x \left (\frac {a}{b}\right )^{\frac {1}{4}} + \sqrt {\frac {a}{b}}\right )}{4 \, {\left (\sqrt {2} a^{2} b^{2} c - \sqrt {2} a^{3} b d\right )}} - \frac {\left (c d^{3}\right )^{\frac {3}{4}} \log \left (x^{2} + \sqrt {2} x \left (\frac {c}{d}\right )^{\frac {1}{4}} + \sqrt {\frac {c}{d}}\right )}{4 \, {\left (\sqrt {2} b c^{3} d - \sqrt {2} a c^{2} d^{2}\right )}} + \frac {\left (c d^{3}\right )^{\frac {3}{4}} \log \left (x^{2} - \sqrt {2} x \left (\frac {c}{d}\right )^{\frac {1}{4}} + \sqrt {\frac {c}{d}}\right )}{4 \, {\left (\sqrt {2} b c^{3} d - \sqrt {2} a c^{2} d^{2}\right )}} - \frac {1}{a c x} \]

input
integrate(1/x^2/(b*x^4+a)/(d*x^4+c),x, algorithm="giac")
 
output
-1/2*(a*b^3)^(3/4)*arctan(1/2*sqrt(2)*(2*x + sqrt(2)*(a/b)^(1/4))/(a/b)^(1 
/4))/(sqrt(2)*a^2*b^2*c - sqrt(2)*a^3*b*d) - 1/2*(a*b^3)^(3/4)*arctan(1/2* 
sqrt(2)*(2*x - sqrt(2)*(a/b)^(1/4))/(a/b)^(1/4))/(sqrt(2)*a^2*b^2*c - sqrt 
(2)*a^3*b*d) + 1/2*(c*d^3)^(3/4)*arctan(1/2*sqrt(2)*(2*x + sqrt(2)*(c/d)^( 
1/4))/(c/d)^(1/4))/(sqrt(2)*b*c^3*d - sqrt(2)*a*c^2*d^2) + 1/2*(c*d^3)^(3/ 
4)*arctan(1/2*sqrt(2)*(2*x - sqrt(2)*(c/d)^(1/4))/(c/d)^(1/4))/(sqrt(2)*b* 
c^3*d - sqrt(2)*a*c^2*d^2) + 1/4*(a*b^3)^(3/4)*log(x^2 + sqrt(2)*x*(a/b)^( 
1/4) + sqrt(a/b))/(sqrt(2)*a^2*b^2*c - sqrt(2)*a^3*b*d) - 1/4*(a*b^3)^(3/4 
)*log(x^2 - sqrt(2)*x*(a/b)^(1/4) + sqrt(a/b))/(sqrt(2)*a^2*b^2*c - sqrt(2 
)*a^3*b*d) - 1/4*(c*d^3)^(3/4)*log(x^2 + sqrt(2)*x*(c/d)^(1/4) + sqrt(c/d) 
)/(sqrt(2)*b*c^3*d - sqrt(2)*a*c^2*d^2) + 1/4*(c*d^3)^(3/4)*log(x^2 - sqrt 
(2)*x*(c/d)^(1/4) + sqrt(c/d))/(sqrt(2)*b*c^3*d - sqrt(2)*a*c^2*d^2) - 1/( 
a*c*x)
 
3.8.84.9 Mupad [B] (verification not implemented)

Time = 10.36 (sec) , antiderivative size = 5962, normalized size of antiderivative = 12.96 \[ \int \frac {1}{x^2 \left (a+b x^4\right ) \left (c+d x^4\right )} \, dx=\text {Too large to display} \]

input
int(1/(x^2*(a + b*x^4)*(c + d*x^4)),x)
 
output
2*atan(((-d^5/(256*b^4*c^9 + 256*a^4*c^5*d^4 - 1024*a^3*b*c^6*d^3 + 1536*a 
^2*b^2*c^7*d^2 - 1024*a*b^3*c^8*d))^(1/4)*(x*(4*a^11*b^9*c^12*d^8 + 4*a^12 
*b^8*c^11*d^9) - (-d^5/(256*b^4*c^9 + 256*a^4*c^5*d^4 - 1024*a^3*b*c^6*d^3 
 + 1536*a^2*b^2*c^7*d^2 - 1024*a*b^3*c^8*d))^(3/4)*(x*(-d^5/(256*b^4*c^9 + 
 256*a^4*c^5*d^4 - 1024*a^3*b*c^6*d^3 + 1536*a^2*b^2*c^7*d^2 - 1024*a*b^3* 
c^8*d))^(1/4)*(1024*a^12*b^12*c^20*d^4 - 4096*a^13*b^11*c^19*d^5 + 6144*a^ 
14*b^10*c^18*d^6 - 4096*a^15*b^9*c^17*d^7 + 2048*a^16*b^8*c^16*d^8 - 4096* 
a^17*b^7*c^15*d^9 + 6144*a^18*b^6*c^14*d^10 - 4096*a^19*b^5*c^13*d^11 + 10 
24*a^20*b^4*c^12*d^12)*1i - 256*a^11*b^12*c^19*d^4 + 768*a^12*b^11*c^18*d^ 
5 - 768*a^13*b^10*c^17*d^6 + 256*a^14*b^9*c^16*d^7 + 256*a^16*b^7*c^14*d^9 
 - 768*a^17*b^6*c^13*d^10 + 768*a^18*b^5*c^12*d^11 - 256*a^19*b^4*c^11*d^1 
2)*1i) + (-d^5/(256*b^4*c^9 + 256*a^4*c^5*d^4 - 1024*a^3*b*c^6*d^3 + 1536* 
a^2*b^2*c^7*d^2 - 1024*a*b^3*c^8*d))^(1/4)*(x*(4*a^11*b^9*c^12*d^8 + 4*a^1 
2*b^8*c^11*d^9) - (-d^5/(256*b^4*c^9 + 256*a^4*c^5*d^4 - 1024*a^3*b*c^6*d^ 
3 + 1536*a^2*b^2*c^7*d^2 - 1024*a*b^3*c^8*d))^(3/4)*(x*(-d^5/(256*b^4*c^9 
+ 256*a^4*c^5*d^4 - 1024*a^3*b*c^6*d^3 + 1536*a^2*b^2*c^7*d^2 - 1024*a*b^3 
*c^8*d))^(1/4)*(1024*a^12*b^12*c^20*d^4 - 4096*a^13*b^11*c^19*d^5 + 6144*a 
^14*b^10*c^18*d^6 - 4096*a^15*b^9*c^17*d^7 + 2048*a^16*b^8*c^16*d^8 - 4096 
*a^17*b^7*c^15*d^9 + 6144*a^18*b^6*c^14*d^10 - 4096*a^19*b^5*c^13*d^11 + 1 
024*a^20*b^4*c^12*d^12)*1i + 256*a^11*b^12*c^19*d^4 - 768*a^12*b^11*c^1...